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We propose a new method for generating new Hamiltonians 
for the development of supersymmetric Quantum Mechanics 
for real as well as complex systems including twins and 
quadruplets Hamiltonians. Interestingly Hamiltonians, 
which can be realised in this method , can hardly be derived 
using all previous formalisms developed so far.Further 
examples discussed here include both shape invariant and 
shape non-invariant potentials .We notice a new thing (which 
has not been revealed earlier) in the case of shape invariant 
potentials i.e by suitable choice of parameter in superpotential 


the generated Hamiltonians satisfy either (i) SUSY conditions 
i.e EM = Et+\ ; Ef> = 0 or (ii) Iso-spectral condition i.e 
E 'jdd = E { n hHowever for non-shape invariant potentials 
energy relations depend only on the choice of 
superpotentials.We present both analytical as well as 

numerical results. 
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I.Introduction 

Generation of Quantum systems has been a subject of main 
interest since the early development of Quantum Mechanics.In 
fact it gained momentum after the work of Witten [1] ,who 
proposed that similar quantum systems can be generated 
through ” Supersymmetry” .However ,Gendeshtein [2] proposed 
that when generated Hamiltonians( H +: HE) reflect similar 
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shape in potentials (commonly known as shape invariant) then 
energy eigenvalues can be related to each other obey only 

SUSY conditions i.e 

4 +) = 43 (i) 

and 

kb = 0 (2) 

where > H + and E^~' —> II. Interestingly .not only a 

report but also a book on supersymmetry has been published 
by Cooper,Klrare and Sukhatme [3] where the authors 
strongly believe in the validity of relations in Eqs(l,2) and 
proposed many exactly solvable Hamiltonians reflecting shape 

invariants potentials i.e 

(x, ao) = V-(x, a\) + R(a\) (3) 

to justify the same.Mean while PT symmetric Quantum 
Mechanics has gained momentum particularly after the work 
of Bender and Boettecher [4] .In order to develop the 
supersymmetry in PT systems Bazeia,Das,Greenwood and 
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Lanso [5] proposed that for shape invariant potentials (as in 
Eq(3)) only Iso-spectral condition 

4 +) = 4 _) (4) 

is the most appropriate one .Of course authors [5] cite 
examples to justify their stand.Here we notice that when the 
generated Hamiltonians are Hermitian in nature they obey 
SUSY condition [3], however when generated Hamiltonians 
determined using siumilarity type of transformations [5] in PT 
symmetric systems ,obey Iso-spectral condition.Hence we feel 
there is lacking co-operation between two groups of authors 
[3,5] .Further from the literature on Supersymmetry we notice 

that Hamiltonian of the type 

H = (p + igf + V(x) (5) 

proposed much earlier by Hatano and Nelson [6] .has not been 
considered [5,3]. Hence from the above discussions we feel that 
existing work on Supersymmetry needs development. Here we 
propose a new method for the development of Supersymmetry 
in Quantum Systems which can be visualized in real or 
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complex Hamiltonians .Futher interesting examples( including 
exactly solvable shape invariant potentials) are presented in 
addition to numerical calculation. 

II.Generators,Superpotentials and Hamiltonians. 

Let us define two generators using operator theory as 

A = ~ + iW(x) (6) 

ax 


B = jL + iW*(-x) (7) 

ax 

where superpotential W(x) can be real or complex.lt may not 
be practically possible to generate either A from B or 
vice-versa. Now using A and B we define two Hamiltonians 

H+ = AB = p 2 -[W(x)-W , (-x)}p-i -- L-w(x)W*{-x) 

ax 

( 8 ) 

and 

H_ = BA = p 2 -\W(x)-W*{-x)\p+i dW } X ^ -W(x)W*{-x) 

ax 

( 9 ) 
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One can notice that AB ^ BA i.e H + ^ H- . This type of 
generated Supersymmetric Hamiltonians are new to the 
literature and can not be visualised by any of the previous 


work[l,2,3,5]. 

IIA.SUSY conditions :Real Hamiltonian 

Let us consider that 

W{x) = ix (10) 

in this case we find 

PTW(x)PT = W(x) = ix (11) 

The generated Hamiltonians are 

H_ = p 2 + x 2 — 1 (12) 

and 

H + = p 2 + x 2 + 1 (13) 

It is easy to see that the corresponding energy levels are 

E ( ~ ] = 2 n (14) 
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( 15 ) 


— 2n + 2 — En +1 
with 

EV = 0 (16) 

where E\P —* //, and E^) —> //_. 

IIB.SUSY conditions :Complex Hamiltonian 


Now consider that 

LH(x) = ix — ig (17) 

In this case for g ^ 0 we find 
PTW(x)PT ^W(x) = ix — ig (18) 

The generated Hamiltonians are 

H_ = (p + fg) 2 + x 2 — 1 (19) 

and 

77+ = (p + ig) 2 + x 2 + 1 (20) 


In order to calculate eigenvalues of the above Hamiltonians, 
we use SU( 11) algebra [6] . Let us consider the Hamiltonian 
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of the type [6] 

H — hup 2 + /i 22 ^ 2 + ihu(xp + px ) + ih\p + h 2 X (21) 


£n— \hnh 22 + hi2(2n + 1 ) + 


having energy eigenvalue 

(/i^/i22 — h\hn — 2h\h2hi2 


4(/in/l22 + hi2 


( 22 ) 


Now using the above formula we find for II ; II { one has 
(i)hi = I 12 = 1 (h) h \2 — h ,2 = 0 (iii) hi = 2 g . Hence it is easy 

to show that 


E( ) — 2 n 


and 


— 2n + 2 — E^+i 


with 


EE = 0 


(23) 


(24) 


(25) 


where E^ —> H + and E { n 1 —y IE. This type of Hamiltonian 
was suggested ealier [7].Interestingly we realise it ina 


7 



systematic way in the context of Supersymmetry and notice 
that SUSY is not broken . Interestingly the above 
Hamiltonian is non -VT invariant and non-Hermiticity in 
nature . Further to support the above example on SUSY 
condition , we consider another interesting form of W ( x ) as 

W (x) = ikx 3 — ix 2 g (26) 

The corresponding Hamiltonians are the following 

H _ = p 2 + 2 igx 2 p — 3 kx 2 + 2 gx — g 2 x A + k 2 x 6 (27) 

and 

H + = p 2 + 2 igx 2 p + 3 kx 2 + 2 gx — g 2 x A + k 2 x 6 (28) 

Now using matrix diagonalisation method [8] we calculate the 
eigenvalues and reflect the same in table-1 Table -1 
First five eigenvalues of SUSY Hamiltonians . 



k 

g 

£(-) 

n 

£(+) 

n 



0 

1.935 482 



1.935 482 

6.298 495 



6.298 495 

11.680 970 



11.680 970 

18.042 635 



18.042 634 

25.254 604 

2 

2 

0 

2.737 184 



2.737 184 

8.907 417 



8.907 417 

16.519 386 



16.519 389 

25.516 139 



25.516 137 

35.715 404 


IIC.Iso-spectral condition :Complex Hamiltonian 

In the previous examples we have considered 
superpotentials having complex in nature.Now we 
consider superpotentials having real in nature i.e 


W(x) = x 2 
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in this case we find 


PTW (x)PT = W(x) = x 2 (30) 

The generated Hamiltonians are 

i7_ = p 2 — x A — 2 ix (31) 

and 

H + — p 2 — x A + 2 ix (32) 

The energy eigenvalues of these Hamiltonians are 
obtained by solving Schroedinger’s equation 
numerically [9] and are tabulated in Table-2.From 
the table one can see that energy levels satisfy the 
iso-spectral condition [7,8] i.e 

4 +) = 4 _) (33) 


Table -2 

First five eigenvalues of Iso-spectral Hamiltonians . 
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Level 

E (+) 

n 

£(-) 

n 

0 

0 

0 

1 

3.398 150 

3.398 150 

2 

8.700 453 

8.700 453 

3 

14.977 808 

14.977 808 

4 

21.999 001 

21.999 001 


IID.SUSY conditions: Exactly solable shape 
invariant Potentials. 

Let us consider one interesting type of complex 

superpotential as 

W (x) = ix — i— (34) 

x 

in this case we find 

PTW (x)PT = W(x) =ix-i- (35) 

X 

The present superpotential contains a linear term 
and an inverse term ,whose similar type was first 
reflected in the work of Supersymmetry by 
Rath[9] in the context of Hermitian 
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Hamiltonian.Later on this type was used by Bazeia 
et.al [5] by making a modification as 


W(x) = ix + e — — 


ix + e 


(36) 


in the context of PT symmetry of supersymmetric 
Quantum Mechanics to discuss iso-spectral 
conditions .Interestingly Bazeia et.al [5] have not 
reflected neither numerical nor analytical results 
to justify their stand. Here the generated 
Hamiltonians are obtained using slight modified 


form of superpotential as 

H+ = P - + X - + A(A + 1 ) -A + 0.5 
+ 2 2 2a; 2 

and 

H_ = — + — + ~ ^ — A — 0.5 

2 2 2a; 2 


(37) 


(38) 


It is easy to see that the corresponding exact 

energy levels are 


E)-) = 2n + 0.5 + 0.5 /l + 4A(A - 1) - A (39) 

E i +> = 2n + 1.5 + 0.5/1 + 4A(A + 1) - A (40) 
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One can see that for any value of A > 2 ,the energy 
levels satisfy the SUSY conditions as 


Ei + \\ > 2) = 2n + 2 (41) 

Et \A > 2) = 2 n (42) 

with 

4 _> (A > 2) = 0 (43) 

Interesting point is that eigenvalues are 
independent of A. The corresponding wave 
functions are the following: 

$(+) = C+[{~ + xf - A(A + l \ n x ^ + ^ + MX + l ) e -^/2 

CLOC OC 

(44) 

and 

= C~[(~ + xf - A(A 7 1) ]"^ 0 ' 5+ V' ft25+A(A - I) e-^ 2 

(45) 


IIE.Iso-spectral condition: Exactly solable shape 

invariant Potentials. 
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Here we consider slight modification of above 
complex superpotential as 


W (x) — ix + i— 
x 

in this case we find 


PTW ( x)PT =W(x) = ix J r i 


(46) 


A 


X 


(47) 


Here the generated Hamiltonians are obtained 
using slight modified form of superpotential as 


p 2 x 2 A(A - 1 ) 

tf + -- + T + ^^ + A + °.5 


(48) 


and 


p 2 x 2 AA+ 1 ) 

II = — H-b —— 77 —~ + A 

2 2 2x 2 
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It is easy to see that the corresponding exact 

energy levels are 


£<-) = 2 n + 0.5 + 0.5y'l + 4A(A + 1 ) + A (50) 

EM = 2 n + 1.5 + 0.5^1 + 4A(A - 1) + A (51) 
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One can see that for any value of A > 2 ,the energy 
levels satisfy the iso-spectral condition as 

£< +) (A > 2) = E ( ~\A > 2) = 2n + 2A + 1 (52) 

Here energy levels are parameter dependent. 

The corresponding wave functions are the 

following: 

= £>+[(_ A + xf - A(A 7 1 ) ]^0-5 + V0^25 + A(A-l) e -^/2 

CLOC OC 

(53) 

and 

$,(-) = D-{(-4- + X ) 2 - ^A±!2] V.5 +Vo^A(ATI) e -^ / 2 
ax x z 

(54) 

III.Generators,Superpotentials : Twins and 
Quadruplets Hamiltonians. 

In the above , we have generated only a pair of 
SUSY or Iso-spectra using suitable choice of 
superpotential . However in this generation we try 
to generate twin SUSY and Quadruplet iso-spectra 
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using the method given below . Let us define two 
generators using operator theory as 

A = —ip — W\ (55) 


B — ip — W 2 



such that the product W\ 115 is an even function of 

x. 

o dWo 

H+ = AB =p 2 + i[W 2 -Wi\p + —P + WiW 2 (57) 

ax 

and 

H-=BA = p 2 + i[W 2 - Wi p - P- + WtW 2 (58) 

ax 

One can notice that AB ^ BA i.e H + ^ i7_ . This 
type of generated Supersymmetric Hamiltonians 
are new to the literature and can not be visualised 
by any of the previous work[l,2,3,5]. Here as per 
the condition stated above one has to select Wp W 2 

even or odd . 
vspace0.2in 
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IIIA. Twins SUSY 


Case-(a) :: W\ = x and W 2 = x 3 

Hi = AB = p 2 + ix 3 p — ixp + 3x 2 + x A (59) 
Hi = BA = p 2 + ix 3 p — ixp — 1 + x 4 (60) 

Case-(b) :: W 2 = x and W\ = x 3 

— AB = p 2 — ix 3 p + ixp + 1 + x A (61) 

Hi = AB = p 2 — ix 3 p + ixp — 3x 2 + x A (62) 

In order to show the twin behaviour , we calculate 
the energy eigenvalues and reflect the same in 
table-3 using the matrix diagonalisation method 

[ 8 ], 

Table - 3 :First four eigenvalues of Twin-SUSY 

Hamiltonians . 
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Ht 

Hr 

m 

H 2 - 

2.067 999 2 

0 

2.067 999 2 

0 

5.631 827 3 

2.067 999 2 

5.631 827 3 

2.067 999 2 

9.995 229 9 

5.631 827 3 

9.995 229 9 

5.631 827 3 

15.047 560 1 

9.995 229 9 

15.047 560 1 

9.995 229 9 


From the table one can see that H and H} have 
the same eigenspectra even though the 
corresponding Hamiltonians are not the same . 
Similarly f/f and H^ pertain to same eigen spectra 
.Hence they are called twin-SUSY. It is worth to 
mention that the generated Hamiltonians are 
complex and satisfy VT - symmetry condition. 

IIIB- Quadruplet Iso-Spectra 
Case-(c) :: W\ = x 2 and W 2 = % 4 

= AB = p 2 + ix A p — ix 2 p + 4x 3 + x 6 (63) 
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and 


H 3 = BA = p 2 + ix^p — ix 2 p — 2x + x 6 (64) 

Case-(d) :: W 2 = x 2 and W\ = x 4 


H 3 = AB = p 2 — ix A p + ix 2 p + 2x + x 6 (65) 

and 

H± — AB = p 2 — ix A p + ix 2 p — Ax 3 + x 6 (66) 

In order to show the quadruplet behaviour , we 
calculate the energy eigenvalues and reflect the 
same in table-4 using the matrix diagonalisation 

method [8]. 

Table - 4 :First four eigenvalues of Quadruplet 
Iso-spectral Hamiltonians . 
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m 

h 3 - 

Hi 

Hi 

0.537 037 9 

0.537 037 9 

0.537 037 9 

0.537 037 9 

4.006 022 7 

4.006 022 7 

4.007 022 7 

4.006 022 7 

9.019 924 8 

9.019 924 8 

9.019 924 8 

9.019 924 8 

15.215 167 0 

15.215 167 0 

15.215 167 0 

15.215 167 0 


In this case the generated Hamiltonians are 


non -VT invariant in nature . 

IV. Conclusion. 

In this paper ,we have developed a new method 
for supersymmetric Quantum System to realise it 
in real or complex operators.For shape invariant 
potentials we show that energy levels are 
intimately related reflecting susy or iso-spectral 
behaviour only on selection of parameters.Hence 
we co-relate previously published two different 
concepts [3,5,10] . We also generate two-Twins 
satisfying SUSY behaviour and quadruplet 
satisfying iso-spectral behaviour .It should be 
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remembered that present selections of A and B are 
not the only possibility . In fact one can follow 
this procedure to generate other types . In order 
to help the reader on new generation , we present 

another type as follows: 

Ai = i-— + W(x) (67) 

ax 

Bx=A + W*(-x) ( 68 ) 

ax 

where superpotential W(x) can be real or 
complex.lt may not be practically possible to 
generate either A\ from B\ or vice-versa.Now using 
A\ and B\ we define two Hamiltonians 

H + = AiBi =p 2 -\W(x)+W*(-x)]p+i - - > -+W(x)W , {-x) 

ax 

(69) 

and 

H_ = B\A\ = p 2 -[W(x)+W , (-x)}p+i dW }K w(x)W*{-x) 

ax 

(70) 
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This new form can easily generate new type of 
SUSY Hamiltonians discussed recently by 
Marqques,Negleni and da Silva [12] . Here we 
select the form of W(x) as 

W ( x ) = i\x\ 2 (71) 

and generate SUSY Hamiltonians as 

H- = p 2 — 2\x\ + x A (72) 

H + = p 2 + 2\x\ + x A (73) 

Further using matrix diagonalisation method , we 
present first five eigenvalues in table-5. 

Table - 5 


First four eigenvalues of New SUSY Hamiltonians 
and previous results [12] . 


£(+) 

n 

UH 

n 

E )+) Previous [12] 

UH 

n 

Previous [12] 

1.969 9 

0 

1.972 3 


0 

5.507 1 

1.969 5 

5.510 0 


1.969 6 

9.394 5 

5.506 8 

9.415 2 


5.508 4 

13.858 3 

9.394 2 

13.893 6 


9.398 6 
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Hence we believe that present development will 
enrich the existing knowledge on supersymmetry 
in the context of energy calculation. 
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